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Sampling and reconstruction of operators 

Gotz E. Pfander, Member, IEEE, and David Walnut 


Abstract —We study the recovery of operators with II 
bandlimited Kohn-Nirenberg symbol from the action of 
such operators on a weighted impulse train, a procedure 
we refer to as operator sampling. Kailath, and later Kozek 
and the authors have shown that operator sampling is 
possible if the symbol of the operator is bandlimited to 
a set with area less than one. In this paper we develop 
explicit reconstruction formulas for operator sampling 
that generalize reconstruction formulas for bandlimited 
functions. We give necessary and sufficient conditions on 
the sampling rate that depend on size and geometry of 
the bandlimiting set. Moreover, we show that under mild 
geometric conditions, classes of operators bandlimited to 
an unknown set of area less than one-half permit sampling 
and reconstruction. A similar result considering unknown 
sets of area less than one was independently achieved by 
Heckel and Boelcskei. 

Operators with bandlimited symbols have been used 
to model doubly dispersive communication channels with 
slowly-time-varying impulse response. The results in this 
paper are rooted in work by Bello and Kailath in the 
1960s. 

Index Terms —Bandlimined Kohn-Nirenberg symbols, 
spreading function, operator Paley-Wiener space, channel 
measurement, channel identification, operator identifica¬ 
tion, operator sampling, Gabor analysis, symplectic ma¬ 
trices. 
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I. Introduction 

In this paper we develop a sampling theory and recon¬ 
struction formulas for operators bandlimited to domains 
of small area. Analogously to the classical sampling 
theory of functions, the objective of operator sampling 
is to fully characterize an object from at first sight 
insufficient information, specifically by observing an 
operator’s action on a single input, typically a discretely 
supported distribution, viz; a weighted delta train. The 
theory developed herein applies to so-called bandlimited 
operators, defined as operators whose Kohn-Nirenberg 
symbol is bandlimited. The symplectic Fourier transform 
of the Kohn-Nirenberg symbol of an operator is referred 
to as its spreading function, so that we are considering 
operators whose spreading function is compactly sup¬ 
ported or is contained in the fundamental domain of a 
lattice. In engineering terms, the operators considered are 
characterized by limited time-frequency dispersion. 


A. Identification and sampling of operators 

The operator identification problem addresses the 
question whether an operator from a given class can be 
recovered from its action on a single probing signal. That 
is, for a given class of operators PL, does there exist an 
input signal g so that Hg determines H. Mathematically 
speaking, we require that the map <f> 9 : H ->• Hg 
be injective on PL. In order to be stable under noise 
introduced, for example, by physical considerations or 
digital processing, it is reasonable to require in addition 
that the map <f> 9 have a bounded inverse |fT6l - 

Definition 1.1: Fet It be a collection of linear opera¬ 
tors mapping a space of functions or distributions X (R) 
to a normed function space Y(R). If for some g £ .A(R), 

P> g : PL —> Y(R), H i—>• Hg 

is bounded above and below, that is, if there are constants 
0 < A < B < oo such that 

A\\H\\ n < \\Hg\\ Y < B \\H\\ n (1) 

for all H £ PL, then we say that PL is identifiable with 
identifier g £ X(R). If PL is not linear, then condition 
dT]) is replaced by 

A\\Hi - H 2 \\ h < \\Hig - H 2 g\\ Y < B ||//i - H 2 \\ H 

( 2 ) 


for all H u H 2 £ H. 

We refer to operator identification as operator sam¬ 
pling when the identifier is a discretely supported distri¬ 
bution. 

Definition 1.2: A strictly increasing sequence A = 
{A n }nez in I is a set of sampling for an operator 
class PL, if for some never-vanishing sequence (d n ) n £Z, 
we have that dn8\ n identifies PL. We define the 

sampling rate of A by 

D( A) = lim 

r—>■ oo T 

where 


n (r) = inf #{n: A n £ [x, x + r]} 

a;SR 

assuming that the limit exists. 

Remark 1.3: (1) D( A) can be interpreted as the aver¬ 
age number of deltas appearing in the identifier per unit 
time and corresponds to the lower Bending density of 
A. 

(2) The assumption that the sequence (d n ) never van¬ 
ishes ensures that the sampling rate depends only on A. 
In particular, we avoid the situation in which for some 
set A' D A, of higher density than A, \' m = 

Jf. n dnd\ r ,. where d' m = d n whenever A' m = A n and 
d’ m = 0 otherwise. 

In this paper we will consider mostly sampling sets 
that are periodic subsets of a fixed lattice on R. 

Definition 1.4: We say that an operator class H can 
be identified by regular operator sampling if there exists 
T > 0, L £ N, and a period-L sequence c = (c n ) such 
that Enez c n d n T identifies H. 

In regular operator sampling, D( A) = |c||o/ (TL) 
where 


||c||o = #{n: 0<n< L— 1 and c n / 0} 

is the support size of the vector (co, ..., cl~ i). In the 
remainder of this paper we will abuse notation and not 
distinguish the vector c £ C L from the doubly-infinite 
/.-periodization c = (c n ). 

Our work addresses the identifiability of classes of 
operators characterized by their Kohn-Nirenberg symbol 
being bandlimited to a set S (the spreading support ). 


B. Operator representations, bandlimited operators, and 
operator Paley-Wiener spaces 

Similarly to linear operators on finite dimensional 
space being represented by matrices, the Schwartz kernel 
theorem implies that linear operators on any of the 
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classical function spaces on R can be represented by 
their kernel, that is, formally, we have 

H fi x ) = J K H (x,y)f(y)dy, (3) 

for a unique kernel K//Q 

As operators are in 1-1 correspondence with their 
kernels, they can also be formally represented by their 
time-varying impulse response h, their Kohn-Nirenberg 
symbol a, or their spreading function //. In fact, formally, 

Hf(x) = J h H {x,t) f(x - t) dt (4) 

= ff VH(t, v) e 2m ^ x ~^ f(x - t ) dv dt (5) 

= J a H (x,Oe 2 ™tmdti, ( 6 ) 

where 


hn(x, t ) = k h(x, x — t) 

= J <xh{x, 0 e 2m ^d^, 

= J rj H (t, v) e 2 ’ Kiv i x ~ t ) dv (7) 

and the Fourier transform in © is normalized as 
Ff{£) = m = j f{x)e- 2 ™tdx. 

Operator representations such as those given in ©, 
©, © are considered in the theory of pseudodifferential 
operators where we write 

cr(x,D)f(x) = j o(xff)e 2nixi f(^)d^ 

With the symplectic Fourier transform given by 

F s F(t, v) = jj F(x, 0 e ~ 2 ^-^) dx d£, 

© implies e~ 2nltv r)H(t,v) = F s aH(t,u). We say that 
the operator H is bandlimited to the set S C R 2 if 
supp t]h = supp F s cjh C S. 

Considering now spaces of such operators we arrive 
at the following definition. 

Definition 1.5: Given a set 5 C R 2 , define the oper¬ 
ator Paley-Wiener space OPW(S) by 

OPW(S) = {H € £(L 2 (R),L 2 (R)): 

supp F s (J H = supp rj H C 5*} 


'in fact, with <S(R d ) denoting the space of Schwartz class functions 
and <S , (R d ) its dual, we can associate to any linear and continuous 
operator mapping S(R d ) to <S / (R rf ) a kernel n £ .S'"(R 2rf ) so that 
0 holds in a weak sense. Below, we shall consider operators acting 
boundedly on the space of square integrable functions L 2 (R) which 
fall in the framework outlined above. We refer to (23 for a more 
detailed functional analytic treatment of operator and function spaces 
involved. 


where £(L 2 (M),L 2 (M)) denotes bounded operators on 
L 2 (R). The space of Hilbert-Schmidt operators in 
OPW(S) is 

opw 2 {S) = opw{s ) n hs{l 2 { k)) 

= {H € £(L 2 (M),L 2 (R)) : 

supp F s (Jh C S, an G L 2 (M 2 )}. 

The reconstruction formulas presented in this paper for 
0PW 2 (S) hold formally for all of OPW(S). Operator 
Paley-Wiener spaces defined by membership of the sym¬ 
bol in generic mixed L p spaces is considered in l24l ; see 
also Section III-DI below for some examples. 


C. Bandwidth of operators and analogies to classical 
sampling of functions 

The terminology operator sampling is intentionally 
suggestive of the classical theory of sampling of ban¬ 
dlimited functions, and is justified for the following 
reasons. 

(1) Classical sampling is in fact a special case of 
operator sampling in the following sense. If for some 
fixed > 0, the operator class 'H consists of operators 
given by multiplication by functions m € L 2 with 
suppm C [—0/2,TT/2], then choosing 0 < T < I/O 
and g = 5 u t, we have that for H e PL corresponding 
to multiplier m, Hg = yf n m(nT)6 n T from which rn 
and hence H can be recovered. In this case, our gen¬ 
eral reconstruction formula (fT8T) reproduces the classical 
reconstruction formula. For details, see Section III-DI 

Finally note that in this case, since ||c||o = L = 1, the 
(operator) sampling rate ||c||o/(TL) = 1/T coincides 
with the sampling rate in the classical sense. 

(2) In analogy with classical sampling, we can give 
a necessary condition on the (operator) sampling rate 
based on a natural measure of the bandwidth of an op¬ 
erator modeling a time-varying communication channel. 
T. Kailath in lfl3Tl defined the bandwidth of a channel 
with spreading function g(t, v) by 

inf{.B > 0: g(t, u) = 0, Vt G R, v ^ (—B/2,B/2)}. 


Taking into account possible gaps in the spreading sup¬ 
port S, we can more precisely define the bandwidth by 


B(S) = sup | supp r](t, v)\ 

tGR 


xs(-,v) do 


( 8 ) 


where xs is the characteristic function of S. This quan¬ 
tity can be interpreted as the maximum vertical extent of 
S. The following theorem gives a necessary condition on 
a set of sampling for the operator class with spreading 
support S. 
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Theorem 1.6: If S is closed and A is a set of sampling 
for OPW 2 (S ) with inf{|A — p\ : A, p, € A} > 0, then 

D( A) > B(S). 

(3) A sufficient condition on the (operator) sampling rate 
is more elusive and is tied to both the area of the spread¬ 
ing support S and its shape. However, if ,Sj is small, then 
it suffices to observe Hg(t ) on a correspondingly small 
portion of the real line. For details, see Section III-EI 

(4) It should be noted that not only is classical sampling 
a special case of operator sampling, but also the well- 
known result that time-invariant operators are character¬ 
ized by their response to a delta centered at the origin. 
Here A = {0} and S’ is a subset of the t-axis with 
B(S ) = 0. For details, see Section HI-D2I 


D. Physical relevance of bandlimited operators 

In communications engineering, <J4]) and ([5]) are com¬ 
monly used as models for linear (time-varying) commu¬ 
nication channels. The time-varying impulse response of 
the channel hjj{x , t) is interpreted as the response of the 
channel at time x to a unit impulse at time x — t, that is, 
originating t time units earlier. Hence, if hn(x,t) f 0 
only for 0 < t < T, then H is causal with maximum 
time-dispersion T. 

If hn(x,t ) = hji{t) then the characteristics of the 
channel are time-invariant and in this case the channel 
is a convolution operator. As mentioned above, such 
channels are identifiable since hn(t) is the response of 
the channel to the input signal 5o(t), the unit-impulse at 
t = 0. 

A mobile communication channel has the property that 
hn(x,t) depends on x, but changes as a function of x 
rather slowly, since the change in the channel, for exam¬ 
ple, by movement of receiver, transmitter, or reflecting 
objects, is slow when compared with the speed of light 
at which information travels. This slow variance can be 
expressed through a bandlimitation of hn{x, t ) as a func¬ 
tion of x, that is, as a support constraint on the spreading 
function of H, = / hn{x,t) dx, as 

a function of u. We conclude that a causal doubly disper¬ 
sive communications channel with maximum time dis¬ 
persion T, and hn(x, t ) bandlimited in x to [—H/2, H/2] 
is represented by a spreading function supported on 
the set [0,T]x [—fl/2, fI/2], that is, by operators in 
OPW 2 ([0, T] x [—fl/2, fl/2]) since \p H \ = \P s <t H \. 

To substantiate this bandlimitation on <jh(x, t) further, 
we denote translation by t by T t : f(x) i->- f(x — t) and 
modulation by u by M v : f(x) i—>■ e 2mux f(xj/T he latter 
is also referred to as frequency shift as M u f = T v f. 


Then © becomes the operator-valued integral 


H= r]H{t, u) T t M u dv dt 


T /-O 


Ph {t, v)T t M v do dt , 

o J-n 

that is, the spreading function is the coefficient vector of 
the time-frequency shifts that a communication channel 
carries out. Hence, OPW 2 ([0, T]x [—fl/2,fl/2]) has 
maximum time-delay T and maximum frequency spread 

n. 


E. Relation to other work 

In 1959, T. Kailath |[T3l . fT4l . Ifl5l asserted that for 
time-variant communication channels to be identifiable 
it is necessary and sufficient that the maximum time- 
delay, a, and Doppler spread, b, satisfy ab < 1 and 
gave a convincing justification for his assertion on signal- 
theoretic grounds. Kailath considers the response of the 
channel to a train of impulses separated by at least a time 
units, so that in this sense the channel is being “sampled” 
by a succession of evenly-spaced impulse responses. The 
condition ab < 1 allows for the recovery of sufficiently 
many samples of hn{x,t) to determine it uniquely. To 
prove necessity, Kailath assumes that the channel can 
be identified by a probing signal, g, essentially both 
time- and band-limited. If the response, Hg, is also so 
limited, the number of degrees of freedom in Hg can be 
estimated. This number is then compared to the number 
of degrees of freedom in the impulse response hn(x,t) 
under the same time and band-limiting assumption as on 
g in each variable. Comparing degrees of freedom leads 
to the necessary inequality ab < 1. 

Kailath’s assertion was given the precise mathematical 
framework described in Section iFAl and proved in fl6l . 

In 1969, P. A. Bello |2| argued that what is important 
for channel identification is not the product ab of the 
maximum time-delay and Doppler shift of the channel 
but the area of the support of the spreading function. 
It is notable that Kailath also asserted something along 
these lines. This means that a time-variant channel whose 
spreading function has essentially arbitrary support is 
identifiable as long as the area of that support is smaller 
than one. Using ideas from fl6l . Bello’s result was 
confirmed in l28l . 

Building on findings in fT6l , ||28l , ||29l a number 
of results have been established that are now part of 
the herein described sampling theory for operators. For 
example, the results in ll28l were extended from the 
setting of Hilbert-Schmidt operators to a much wider 
class of pseudodifferential operators in lf24l . In fl2l . the 
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choice of non periodic (irregular/jitter) sampling loca¬ 
tions for operator sampling was discussed. Necessary and 
sufficient conditions for the identifiability of bandlimited 
Multiple Input Multiple Output (MIMO) channels were 
given in ||2§1l . 

More recently, sampling results for stochastic oper¬ 
ators, that is, for operators with stochastic spreading 
functions, have been obtained lf22ll . lf32ll . OH . Also, in 
applications, it is required to replace the identifier con¬ 
sidered in this paper by finite time or finite bandwidth, 
that is, smooth, signals. Local recovery results in this 
setting, as well as a reconstruction formula that allows 
for the application of coarse quantization methods prior 
to the approximate recovery of the operator are given in 
Ifl8t Focusing on a parametric setup, the identification 
of bandlimited operators was analyzed with respect to 
applicability in super-resolution radar (TJ. 

In Section HI-FI we address the problem of iden¬ 
tifiability of operators with unknown bandlimitation. 
Independently, Fleckel and Boelcskei ( 0 , 0 ) have 
obtained a result similar to Theorem 12.251 characterizing 
identifiability of a certain (nonlinear) class of operators 
with spreading support of area < 1/2. Theorem 12.251 
gives a sufficient condition for a more general class of 
operators, and Theorem l2.26l generalizes the result in 0 , 
0 . In addition, Fleckel and Boelcskei ( 0 , 0 ) prove a 
remarkable result in which they prove identifiability for 
unknown support sets of area less than one, rather than 
< 1 / 2 . 

II. Main results 

A. Properties of Gabor system matrices 

The basic strategy for operator sampling described in 
this paper was laid out in lf28l . The idea is to translate the 
reconstruction problem into an a priori under-determined 
linear system whose coefficients come from a finite 
Gabor system, and then give conditions under which 
that system can be solved. More specifically, given 
H € OPW 2 (S), T > 0, and L € N, let g = c n S nT 
for some period-L sequence c = (c n ). Then from the 
response Hg(x), we can derive the Lx L 2 linear system 

Z Hg(t, v) = G(c) r/ H (t , u) (9) 

where Z# 9 (t, v) is an /.-vector computed directly from 
Hg, r] H (t,u) is an L 2 -vector consisting of shifts of a 
periodized version of the spreading function rju of H 
(see Lemma 13.71) . and G(c) is an L x L 2 Gabor system 
matrix defined as follows. 

Definition 2.1: Given L € N, let u = e 2m / L and de¬ 
fine the translation operator T on (xo, ..., xl-i) € C L 
by 

Tx = (xl_i,x 0 , Xi, . . . ,XL- 2 ), 


and the modulation operator M on C L by 

Mx = (w°Xo, CJ 1 ®!, . . . , 0 J L ~ 1 XL-i). 

Given a vector c € C L the finite Gabor system with 
window c is the collection {T q M p c} q “1 0 . Define the 
full Gabor system matrix G(c) to be the Lx L 2 matrix 

G(c) = [ D 0 W L | D 1 W L | • • • | D l _ x W l ) (10) 

where D *. is the diagonal matrix with diagonal T k c = 
(■ CL-k , ■ ■ ■, cl- i, c 0 , ..., c L - k - 1 ), and W L is the LxL 
Fourier matrix Wl = (e 2 ™ m / i )^) m 1 =0 . 

Remark 2.2: (1) For 0 < q, p < L — 1, the (q + 
l)st column of the submatrix D p Wl is the vector 
M p T q c where the operators M and T are as in Defi¬ 
nition GLU This means that each column of the matrix 
G(c) is a unimodular constant multiple of an element 
of the finite Gabor system with window c, namely 
{ e - 2 nipq/L T q M p c} L q ~l 0 . 

(2) Note that the finite Gabor system defined above 
consists of L 2 vectors in C L which form an overcomplete 
tight frame for C L iGOl . For details on Gabor frames in 
finite dimensions, see |[20l . |fT9j , |9| and the overview 
article ll33l . 

The reconstruction formulas in this paper are based on 
explicitly and uniquely solving ©. for this purpose we 
require conditions on G(c) under which this is possible. 

Definition 2.3: (§1 The Spark of an M x N matrix F 
is the size of the smallest linearly dependent subset of 
columns, i.e., 

Spark(F) = min{||x||o: Fx = 0, x 0} 

If Spark(F) = M + 1, then F is said to have full Spark. 
Spark(F) = k implies that any collection of fewer than 
k columns of F is linearly independent. 

The existence of Gabor matrices with full Spark has 
been addressed in |[20l and |[2Til . 

Theorem 2.4: lf20l If L is prime, then there exists a 
dense, open subset of c E C L such that every minor of 
G(c) is nonzero. In particular, for such c, G(c) has full 
Spark. 

Note that if L is not prime then the result of this 
theorem does not hold. That is, if L is not prime, then 
for any c € C L there is a minor of G(c) that vanishes. 
However, it has recently been shown by Malikiosis that 
for any L € N, we can get the second half of the 
conclusion. 

Theorem 2.5: lf2TB For every L € N there exists a 
dense, open subset of c € C L such that G(c) has full 
Spark. 

This next result states that, again assuming that L is 
prime, the Spark of the matrix G(c) is related to the 
support size of the vector c. 
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Theorem 2.6: |30l If L € N is prime, and k < 
L, there exists c € C L with the property that 
Spark(G(c )) = k+1, and supp(c) C {0, 1, ..., k— 1}. 
Moreover, the set of such c forms an open, dense subset 
of C k x {0}. 

These theorems show that it is possible to choose a 
period-L sequence c such that the system © always has 
a solution as long as there are no more than L non¬ 
vanishing unknowns on the right side. In fact, if L is 
prime, we can say a bit more, namely that if there are no 
more than k < L non-vanishing unknowns on the right 
side, then we can guarantee solvability with a window c 
supported on no more than k contiguous indices. 

B. Necessary and sufficient conditions for identifiability 
of OPW 2 {S) 

In this section, we explore conditions under which 
the operator class OPW 2 (S) is identifiable. We give 
necessary and sufficient conditions on S under which 
identification is possible with any identifier, then char¬ 
acterize when this is possible using regular operator 
sampling. 

In lfl6l . l28l (cf. lf29l and ||24|), the following result 
is given. Here and in the following, IS 1 ! denotes the 
Lebesgue measure of the set S. 

Theorem 2.7: OPW 2 (S) is identifiable by regular 
operator sampling if S is compact and ,3j < 1, and 
not identifiable if S is open and |Sj > 1. 

The following result guarantees the existence of a 
discretely supported identifier for support sets S with 
| S' | < 1 that satisfy certain periodization conditions. The 
result characterizes operator Paley Wiener spaces that 
can be identified by regular operator sampling. 

Theorem 2.8: Let g = c n d n T with c € C L 

chosen so that G(c) has full Spark. For SCR 2 the 
following are equivalent. 

(i) The map <& g : OPW 2 (S) -> L 2 (R), H Hg is 
injective. 

(ii) The function g identifies OPW 2 {S). 

(iii) S is a subset of a fundamental domain of the lattice 

LT7L x (1/T)Z, that is, 

Xs+(kLT,e/T) < 1 CL.e. (11) 

k,t 

and S periodized by the lattice T7L x 1/(TX)Z is 
at most an L- cover, that is 

Xs+(kT,e/(TL)) < L a.e. (12) 

k,t 

See Figures [TTl3l for an illustration of spreading sup¬ 
ports sets S that lead to identifiable operator Paley 
Wiener spaces. 


v 



T LT t 


Fig. 1. The space OPW 2 (S) is identifiable for S (in blue) with area 
1 as it clearly satisfies CCD and {©■ S has a ( T , 3)-rectification and 
B(S) = 1 i. Such sets were considered in (J_6j, HD, ii). Recovery of 
operators in OPW 2 (S) is possible using the reconstruction formula 

Gl. 


V 



Fig. 2. The union of the colored sets, S, satisfies GD and O- 
Flence, OPW 2 (S) is identifiable by a weighted delta train with 
period-3 weighting sequence even though a ( T , 3)-rectification is not 
possible (note that 7 > 3 boxes are active). Recovering g from Hg 
using ((© directly requires solving three systems of linear equations, 
one to recover g on the yellow support set, one to recover g on 
the red support set, and one to recover g on the blue support set. 
H £ OPW 2 (S) can be reconstructed using formula 1 1 71 . Note also 
that B(S) = 212 and that the sampling rate is 1/T = 312 > 212. 

Remark 2.9: (1) It is clear that if S is bounded, then 
(fill is satisfied as soon as S is contained in a rectangle 
of width TL and height 1/T. 

(2) Note that (fl2l) implies that |Sj < 1, and that if |Sj = 
1, the cover must be an exact T-cover, that is, 

XS+(kT,£/(TL)) = L a - e - 

k,e 

(3) As discussed in detail in Remark 12.121 below, for any 
compact set S with |Sj < 1, there exists T, L so that 
(fTH and (fT2l) hold. 

(4) Note that (fill) and (fl2l) are satisfied for some 
unbounded sets with area less than or equal to one, for 
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example, 

OO 

5= ( U [rz + 1 - 2-^, 7Z + 1 - 2-^ +1 >]) x 

n=0 

permits the choice of T = 1 and L = 1. 

(5) On the other hand, it is not hard to construct an 
unbounded set S of arbitrarily small measure so that 
for all T and L, (fill) fails. Indeed, let { q n }neN be an 
enumeration of the countable set of rational numbers Q. 
For e, <5 > 0 set 

OO 

e,<5 = [ e,e] LJ [_J 2 n [—e,e\+q n S j x [—(5,(5]. 

n=0 

We have < 8e<5 since we are taking the union of 
sets that are not disjoint, in fact, every set in the union 
contains countably many sets in the union. In order to 
show that there exist no T > 0 and L € N such that 
CD holds, observe first that clearly, LT > e, and there 
exists no € N so that \q no — LT\ < e. But then S et s~LT 
intersects with S e> $ on a set of positive measure since 

| [-e, e] n 2 -n [—e, e] + q no - LT\ > 0. 

(6) If Spark{G{c)) = K < L, then OPW 2 (S ) is 

identifiable if the upper bound L in (fl2l) is replaced by 
K — 1. 

(7) The conditions (fTTb and (1T2T) are related to the 
rectification of the region S, that is, its efficient covering 
by small rectangles. (See Fig [3]). 

Definition 2.10: Let S C M 2 , |Sj < 1, T > 0, 
and L e N be given. We say that S admits a (T, L )- 
rectification if 

(a) S is contained in a fundamental domain of the 
lattice (TL) Z x (1/T)Z, and 

(b) the set 

S° = [J S + (kTL } £/T) (13) 

(M)ez 2 

meets at most L rectangles of the form R q , rn = 
[0, T] x [0,1/TX] + (qT, m/TL), 0 < q,m < L. 
The active boxes in the rectification are indexed 
by 

T = {(q,m), 0 <q,m<L: R q)m nS° 0}. 

It is clear that (fTTb and (ITTb are satisfied if S admits 
a (T, L)-rectification, but Figure [2] illustrates that the 
converse is not true. Flowever, {ED and (flTb allow for 
the lineal' system © to change depending on the point 
(f, a). In fact, such an observation further characterizes 
regions S such that OPW 2 (S) can be identified by 
regular operator sampling. 


Lemma 2.11: Suppose that for some T > 0 and L € 
N, S satisfies (flD . Then 

y! Xs+(kT,£/(TL)) < K < L a.e. 
k,e 

if and only if there exists a partition {Aj}^ =l of the set 
[0, T ] x [0,1/(LT)] with the property that for each j, at 
most K of the sets Aj + (kT,£/(LT)), 0 < k,£ < L 
meet S°. Moreover, S can be partitioned as 

N 

s = U < 14 > 

i=i 

where 


S, 


sn 


U Aj + (kT, £/{LT)) 

- k,l£Z 


and where each Sj admits a (T, L)-rectification with 

|r| < k. 

Remark 2.12: (1) If S C M 2 is compact and | S'! < 
1, then it is always possible to choose T > 0 small 
enough and L S N large enough that S admits a (T, L)- 
rectification. In fact we can also require that for all e > 0 
sufficiently small, 

= 5Z 1-^9, m\ < 1-51(1 + e) < 1. 

(q,m)G r 

(2) Under certain mild regularity assumptions on a 
domain S, we can explicitly estimate T and L that work. 
Specifically, L E N can be chosen so that all such 
domains have a (y/L, L)-rectification. 

Theorem 2.13: Fix A, B, e, U > 0, N € N, 0 < a < 
1. Suppose that S C [— A, A] x [— B, B\ and there exist 
N Jordan curves C t such that 

1) S is contained in the interior sets of the Jordan 
curves, 

2) the sum of areas of the interior sets is less than 
a — e, and 

3) the sum of lengths of the Jordan curves is bounded 
by U. 

Then for every L satisfying A, B < (L — l)/2 and 
4:(U/VL + N/L) < e, the set S+ {A, B) has a (y/L, L)- 
rectification with ITI < aL. 


C. Sampling and reconstructing operators 

One of the contributions of this paper is to give explicit 
reconstruction formulas for the impulse response of the 
channel operator from the operator’s response to the 
identifier. Such formulas illustrate a connection between 
operator identification and classical sampling theory. 
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Fig. 3. The set S in blue, its rectification in gray. We have L = 9 
and TO = 1/9. 


1) Operators with rectangular spreading domains: 
We begin by recalling a result from Il24l . It is a special 
case of Theorem 12 .15 1 below, and is the simplest example 
on how Shannon’s sampling theorem can be extended to 
apply to operators. 

Theorem 2.14: For H € OPW 2 {S), S C 
[0, T)x [—fl/2, fl/2) compact and TO, < 1, 


h(x, t ) = e-* u / T J2 {H E M (* + nT ) 
n£Z ^ fc£Z 

sin(^((x — t) — nT)) 


7 r((x — t) — nT) 


*[o,T](f). (15) 


where the sum converges in L 2 (M 2 ) and for each t, 
uniformly in x. 

2) Non-rectangular, rectifiable spreading domains: 
The following theorem gives a reconstruction formula 
for operators in OPW 2 {S) when S has a rectification 
in the sense of Defintion 12.101 

Theorem 2.15: Suppose that S C M 2 and that for 
some {t 0 , n 0 ), S — {to, no) admits a (T, L)-rectihcation, 
and let 0 = 1/{TL). Then OPW 2 {S) can be identified 
by regular operator sampling, and there exist period-L 
sequences b^ m) = and functions $( 9 , m ) (f, z') 

for {q, m) € T, such that 


h{x,t) = e 2 ^ t+to >° 

E E ihq,m),k Hg{t - {q- k)T) 

k 

e — 2 irim(q—k)IL ^ (*, x _ (t + t()) + (g _ k)T) ] . 

(16) 


where the sum converges unconditionally in L 2 (M 2 ). 
Here 

(t, s) = J e 2mus X S(#im) (■ t, n) dn 


where 

S( q , m) = Sn (J {Rq,m + {k/0,£/T)). 

k,lG Z 

Remark 2.16: (1) The coefficient sequences are 
defined in (PTI) and are the rows of a left-inverse of the 
L x |T| submatrix of G{c) that allows i© to be uniquely 
solvable, extended to have period L. 

(2) In light of Lemma 12.1 11 it follows that for any 
region SCI 2 for which regular operator sampling of 
OPW 2 {S) is possible, a formula like (flhl) holds. By 
realizing S' as a disjoint union of sets Sj as in ( 1 1 41 ). each 
of which admits a (T, Lj-rectification. and moreover 
where each {t,n) € Sj corresponds to the same reduced 
linear system in we can write 

N N 

hit: hit, v)xs d {t, = hj{t, v) 

1 =i i=i 

and by © 



For each j, we can take to = zd = 0 in (fl6l) and obtain 


hj {X, t)=J2 J2 Hq,m),k H Si t ~ (9 - k)T) 
k ( q,m)GV 0 

e —2nim(q—k) IL ^ (u _ t+(? _ fc)T)] 

where F ? indexes the active boxes in the (T, Lj- 
rectification of Sj, 

and 

S Lm)= S i n U {Rq,m + {k/0,£/T)). 

k/G Z 

Setting = 0 if {q,m) £ T jt 

N 

h{x,t) = hj{x, t) 

3 =i 

N L -1 

= EE E H q ,m),k H sit-iq-k) T ) 

j= 1 k q,m= 0 

e -2mm{q-k)/L ^^ x _ t + {g _ fc)T) ] 

L—1 

= EE Hg{t-{q-k)T) 

k q,m= 0 

®(q,m),k{t, x — t + {q — k)T) (17) 

where 

N 

3 = 1 






























3) Smooth reconstruction functions in the “oversam¬ 

pled” case: Note that Theorem 12.141 and Theorem 12.151 
both involve the use of sharp cut-off functions in the 
definition of the reconstruction functions s). 

The passage to smooth cut-off and hence reconstruction 
functions is enabled by the assumption that S is compact 
with |5| < 1. This allows for faster decay of the recon¬ 
struction functions, and for the validity and convergence 
of the reconstruction sums in more general function 
spaces. These matters have been studied extensively in 
ll24l . Specifically, we have the following generalization 
of Theorem 12.151 

Theorem 2.17: Suppose that S C R 2 , |5| < 1, is 
compact. Then there exist T > 0, L G N, (to, vf), and 
a period-L sequence c = (c n ) such that g = '/T jn c n 5 n T 
identifies OPW 2 (S). Moreover, there exist period-L 
sequences b^ qrn p (q,m) G T such that 

h{x, t ) = e 2 ^ t+to)u ° 

E E [ b (q,m),k Hg(t -{q- k)T ) 

k (q,m)G r 

e 2 nim(x-t)/LT ^ + * Q ) + ( g _ k ) T ) r ( t _ qT )\ 

(18) 

where r, 4> G S (R) satisfy 

^r(t + fcT) = l = ^( 7 + n/LT), (19) 

fcez nGZ 

where r(t)fi{ 7 ) is supported in a neighborhood of 
[0,T]x[0, 1/LT], and where the sum in (fj~ 8 l) converges 
unconditionally in L 2 and for each t uniformly in x. 

Equation (fl 8 T) is a direct generalization of (fl5T) under 
the assumption that r{t) = Xr 0 i t](^) and ( 2 ( 7 ) = 
x [o,a] (7)- 

4) Rectification by parallelograms: It can be advan¬ 

tageous to consider S' to be a subset of a fundamental 
domain of a general lattice All 2 where A = ( all ) ■ 
Our next theorem relies on basic insights on the role of 
symplectic geometry in time-frequency and generalizes 
Theorem 12.151 For simplicity, we restrict our attention 
to lower triangular matrices (““ ). In Section IIII-II 

we discuss the general case in detail and compute the 
quite involved resulting reconstruction formulas (l38l> - 

<ED. dH, d46ji. 

Theorem 2.18: Let S C R 2 , |S| < 1, and assume that 
with A = (^ ^), det A = Tfl = l/L, for some vq € R, 
S + ( 0 , vq) is contained in a fundamental domain of the 
lattice LAI 2 , and that ^ 0 )) admits 

a (T, L)-rectification, that is, if P q , m = ^4.([ 0 , 1] 2 + 


(q, m) T ), q , m G Z, then 

M([ 0 ,lf)n[JS+( 07 o)+M(M) T C IJ P q ,m. 

( 2°) 

Then OPW 2 (S) can be identified by operator sampling. 
Namely, with the period-L sequences c = (c n ) and 
from Theorem 12.151 and functions 

$( 9 ,m) (t, s) = j e 2mi/s X S(q _ m) (t, v) do, 

S( q ,m) = S n U (P q , m + LA(k,£) T ), 

k/GZ 

h(x,t)=e~™ t2 / T J2 E K,rn),k^ iant,T - {q - k)) 

k ( q ,m)G T 

Hg(t ~(q- k)T) x-(q- k)T) 

( 21 ) 

Here the identifier g = c n e mTan2 5 n T and the recon¬ 
struction sum converges unconditionally in L 2 (R 2 ). If 
the product Ta is rational, say Ta/2 = p/q in lowest 
terms, then ( c n e mTan2 ) n is periodic with period being 
the least common multiple of q and L. In particular, if 
LTa/2 is an integer, then the period is L as well. 

Example 2.19: (1) Figure 0] illustrates Theorem 12. 181 
In this case, S is the union of the red and yellow 
triangles and hence is a parallelogram of area 1 and 
A = (^q) with Tfl = l/L = 1/3. Theorem 12.181 
says that OPW 2 (S) can be identified by a periodically 
weighted delta train of period 2 L = 6. However, since 

{ T o l)A~ 1 S={_\° 1 )S = [0,LT]x[0,n\ 

admits a (T, L)-rectification with L = 3, recovery of the 
spreading function would only require solving a single 
3x3 linear system or, equivalently, finding the three 
period-3 sequences 6 ( g>m ) in (I2T1) would require inverting 
a single 3x3 matrix. 

(2) Alternatively, by considering the red and yellow 
regions separately as in Remark 12.161 21. OPW 2 {S) 
can be identified by a periodically weighted delta train 
of period 3. However, recovery of rj(t, v) requires the 
solution of two 3x3 1 i near systems and finding the 
coefficients m . in (fT71) requires inverting two 3x3 
matrices. 

Example 2.20: Figures [5] and [6]illustrate a situation in 
which OPW 2 (S) can be identified by operator sampling 
but not by regular operator sampling. In this case, 
^ = ( 72 ^ 2 + 1 / 2 )’ = an< ^ hence = 

[0, l] 2 admits a (T, L)-rectification with T = L = 1. 
Therefore, following the notation in the proof of The¬ 
orem 12.181 A = B, c n = 1 for all n, L' = 2, and 






















10 



Fig. 4. The space OPW 2 (S) where S is the union of the red 
and the yellow sets is identifiable with reconstruction formula CD 
with a periodically weighted delta train of period 6. Alternatively, we 
can identify OPW 2 (S) using a periodically weighted delta train of 
period 3, but we have to solve 2 linear systems or equivalently invert 
2 3x3 matrices. In this case reconstruction is given by CD. See 
Example 12. 191 for details. 



Fig. 5. (a) The the operator class OPW 2 (S) with S = 

(2, 2 ; v/2, x/2 + 1/2)[0, l] 2 whose area equals 1 and bandwidth 
equals 1/2 is identifiable by a (non-periodically) weighted delta 
train with sampling density 1/2. It is not identifiable using regular 
operator sampling. (b) T = 1 periodization of S. For details, see 
Example 12.201 


c' n = 1 — e mn . By equation (1441) . OPW 2 (S) can be 
identified by 

g = p{B)g = - e ™)e™ 2 ^ 2 5 n , 

a delta train with non-periodic weights. 

Note that B(S) = 1/2, and that since 1 — e mn = 0 
when n is even, the sampling density of the identifier 
g is also 1/2. Therefore, by Theorem 1 1.61 this identifier 
achieves the minimal sampling rate for this region. 

Next we observe that this region cannot be identified 
by regular operator sampling for any value of T or 
L. Since |S'! = 1, by Remark 12.91 2). the T7L x f YL- 
periodization of S must be an exact L- cover. In other 
words, the inequality in (fl2l) must be an equality. It can 
be shown, however, that for any value of T and L , this 
is not possible. Details of the argument can be found in 
Section IIII-JI 


r 

A 


t 1.85.. 2.12.. 

(b) 

Fig. 6. (a) For periodic operator sampling to succeed with S having 

area 1, we require that the T, ri periodization of S leads to an exact 
L cover of the time-frequency plane. (b) The central piece of the 
set S. For the significance of this set, see Example 12.201 


V 



(a) 


spaces OPW°°’ 2 (S) and OPW 2 ’°°(S) by 


D. Operator Sampling as a Generalization of Classical 
Sampling. 


OPW°°' 2 (S) = {H € C(L 2 (R),L 2 ( K)): 

supp p H C S, \\ct h \\l°°^ < oo} 


By generalizing the setting to other function spaces, 
we can more precisely illustrate the connection between 
operator sampling and the classical sampling theorem 
of Shannon, Whittaker, and Kotelnikov among others, 
and also the connection with the well-known fact that a 
time-invariant operator can be identified by its impulse 
response. 

Definition 2.21: We define the operator Paley-Wiener 


where 


\&h\\l° 


and 


0Tt(t£) \ 2 d£ 


1/2 

OO 


OPW 2 ’°°(S) = {H € £(L 2 (M),L 2 (M)): 

supp ph O S, ||cr//||zA~ < oo} 
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where 

2 \ 1/2 

O’tffoOHotAj 

g|24ll . Theorem 4.2). OPW p,q (S ) is a Banach space with 
respect to the norm \\HWopw^^ = 

Note that convolution with a compactly supported 
kernel whose Fourier transform is in L 2 is an operator in 
OPW 00 ’ 2 and multiplication by a bandlimited function 
in L 2 is an operator in OPW 2, °°. 

1) Identification of convolution operators: First, take 
H to be ordinary convolution by hn(t), that is, 
hn(x,t) = hn{t). In this case H can be identified in 
principle by g = 6 q, the unit impulse at the origin, since 
Hg(x ) = hn(x). That is, A = {0} is a sampling set for 
the class of convolution operators. Translating this into 
our operator sampling formalism results in something 
slightly different. 

Assume that h is supported in the interval [0, T'], 
h £ L 2 and that T > T', and 12 > 0 are chosen 
so that 12 T < 1. In this case, r]H(t,u) = h(t) So(u) 
and (Jh(x,C) = Therefore an £ L°°’ 2 and 

H £ OPW°°' 2 {l 0, T '}x [-12/2,12/2]). 

Applying Theorem 12.171 to this situation, note that 
if 9 = S n T then Hg is simply the T-periodized 
impulse response h(t), and it follows from the theo¬ 
rem (or by direct calculation) that with r, <p £ 5(K), 
r(t) = 1 on [0, T'] and vanishing outside an interval 
of length T containing [0, T'], and with (p( 0) = 1 and 
supp^C [—12/2,12/2], 

+ kT)p(x — t — kT) 
fcez 

= ££ r(t) h(t + kT — nT) ip(x — t — kT) 

fcSZnSZ 

= ^2 h(t) <p(x — t — kT) = h(t). 
fcez 

Flere we have used the fact that r(t) = 1 on [0, T'] 
and vanishes outside a neighborhood of [0, T'] and that 
ffk t( x — t — kT) = 1 by the Poisson Summation 
Formula and in consideration of the support constraints 
on (p. Indeed the theorem says that the sum ip(x — 
t — kT) converges to 1 in the L°° norm and in particular 
uniformly on compact sets. 

2) Identification of multiplication operators (Classical 

Sampling): To compare Theorem 12.171 with the classical 
sampling theorem, take H to be multiplication by some 
fixed function m € L 2 with suppm C [—12/2,12/2] 
then (t)m(u), h(t,x) = 5o(t)m(x — t), 

and an(x,^) = m(x). Let 12' > 12 and T > 0 be 
such that 12'T < 1. Then an € L 2,0 ° and H £ 
OPW 2 '°°([—T /2, T /2] x [-12/2,12/2]). 



Choose r, p £ 5(M) such that suppr C [—T/2,T/2], 
r(0) = 1, supp^ C [— 12'/2,12'/2], and p(u) = 
1 on [-12/2,12/2], If g = then Hg = 

m(nT) 5 n T, and it follows from Theorem 12.171 (and 
by direct calculation) that 

So(t) m(x — t) 

= >-«)£ (Hg)(t + kT)p(x — t — kT) 
fee z 

= r(t) ^2 ^2 m ( nT ) - t-kT) 

fcSZ nSZ 

= ^2 minT) <p(x — nT) 

nGZ 

by support considerations on the function r(t). Therefore 
we have the summation formula 

m(x) = *22 tn(nT) p(x — nT) 
n^L 

where the sum converges unconditionally in L 2 . This 
recovers the classical sampling formula when sampling 
above the Nyquist rate. 


E. Sufficient conditions on the sampling rate in operator 
sampling 

As was observed earlier, a natural measure of the 
sampling rate in operator sampling is the quantity 72(A) 
(Definition 11.21 ). which in the case of regular operator 
sampling is ||c||o/(2~X). A necessary condition on the 
sampling rate in operator sampling was give in terms 
of the bandwidth of a channel (Theorem 11.61 ). The goal 
of this subsection is to investigate sufficient conditions 
on the sampling rate in regular operator sampling that 
guarantee identifiability. 

In the classical sampling theory of functions, the 
sampling rate must exceed the reciprocal of the area of 
the bandlimiting set; and regardless of the measure of 
the bandlimiting set, a (possibly high density) sampling 
set always exists. As mentioned above (Theorem 12.71) . 
operator sampling of OPW 2 (S) is only possible if the 
measure of S satisfies | .S'] < 1 , and necessary sampling 
rates in operator sampling depend on the geometry of S. 

The main result in this paper relevant to finding a suf¬ 
ficient condition on the sampling rate for identification 
of OPW 2 (S) is the following. 

Theorem 2.22: Let S C R 2 be compact, |5| < 1, 
e > 0, and suppose that S has a (T, iV)-rectification 
satisfying (|T| + 2)/A^ < |5|(1 + e) < 1. Then for 
every sufficiently large L £ N, OPW 2 (S) can be 
identified via regular operator sampling by an identifier 
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g = Yh n c n 5 n T, where c = (c n ) is a period-L sequence 
satisfying 



Moreover, if L is prime, then c can be chosen with Cj = 0 
if ||c||o < j < L, that is, such that c is supported on its 
first ||c||o indices. 

Remark 2.23: (1) Note that once an appropriate 
(T, N) -rectification of S is found, the parameter T 
associated to that rectification is fixed. Subsequently, 
a periodic weighting sequence can be found for the 
delta train Cn^nT whose relative support is bounded 
essentially by the area of S. Moreover, if L is prime, 
and c is supported on {0, 1, ..., ||c||o — 1}, then this 
represents a bunched operator sampling that can allow 
for the efficient identification of the channel in the 
following way. 

If the area of S is small, and if K represents the 
“memory” of the channel (that is, for each v, ij(t. v) 
is supported in the interval [0, K ]), then the response of 
the channel to the delta train c nfnT is supported on 
the set 

U {[0,T\\c\\o + K\+jLT) 
je z 

and hence vanishes on the set 

U {[T\\c\\o + K,LT\+jLT). 

jez 

The “dead time” represented by this set can be used for 
other purposes. Note also that | LT — (T||c||o + K)\ > 
LT( 1 — |S|(1 + e) — K/(LT )) so that the length of the 
dead time within each period of the channel response 
increases with L. 

(2) Another interpretation of this result is that the sparsity 
of the matrix G(c) in the linear system © can be 
controlled by the area of the spreading support. In this 
case, ||c||o/T gives the fraction of nonvanishing entries 
in each column of G(c). Hence S with small support 
guarantees that G(c) can be chosen to be sparse. 

Remark 2.24: The “dead time” referred to above can 
be thought of as a measure of the capacity of the 
unknown, bandlimited channel in the sense that only 
during this time can data be sent over the channel. With 
this notion of capacity, the above discussion says that 
the capacity of a time-varying channel decays linearly 
with the area of its spreading support. 

F. Sampling and reconstruction of operators with small, 
but unknown support 

Just as in classical sampling, operator sampling re¬ 
quires full knowledge of the bandlimitation we expect an 


operator to have, that is, the reconstruction formulas for 
OPW 2 (S) depend on knowing the region S. However, 
in some applications S may not be known precisely, but 
only some information on its size, geometry and location 
is given by physical considerations. In this section we 
address the question whether such operators can be 
sampled and reconstructed in a stable matter. 

Theorem 2.25: For A,B,U,e,cr > 0 and N € N, 
let 'H(A, B,U, N, e,a) contain all operators such that 
supp F s cth = supp ph C [—A,A]x[—B,B\ satisfies 
the hypothesis of Theorem 12.131 with a < 1/2. Then 
there exists L e N and an /.-periodic sequence (c n ) such 
that g = J2n c nS n /fL identifies W.{A, B, U, N, e, a). 

The reconstruction of an operator H € 
H(A, B,U,N,e,cr) is then carried out as follows. 
First choose L as in Theorem 12.131 and let Rh denote 
the rectified support of H, that is, the union of 
(l/VX) x (1/vT) boxes that cover supp tjh having 
area not greater than 1/2. Under this assumption, we 
determine Rh- In the final step, we apply the operator 
reconstruction formula developed in Theorem 12.151 to 
OPW 2 (R h ). 

To determine the rectified support of tjh with H € 
H(A, B,U,N,e,cr), we will apply ideas from com¬ 
pressed sensing. Indeed, Lemma 13.71 below, shows that 
from H Yh n C n& ny /Z’ we can compute a length L vector 
y(t, v) with y (t, v) = G(c)x(f, v) and where the un¬ 
known discrete support of the length L 2 vector x(t, v) 
encodes the support of the bivariate function r)jj(t,v). 
In fact, recovering the vector x(t, v) for a single point 
(t, v) provides us with the support structure of r///. Note 
that the conditions given above imply that x(t, v) has at 
most L/ 2 nonzero components. 

The full-Spark matrix G(c) plays the role of a mea¬ 
surement matrix and has the ability to recover any L/ 2- 
sparse vector x(f, v) lfl9l . ll20ll . But finding an L/2- 
sparse vector requires consideration of every support 
structure out of (l/2) possible ones, which is hardly 
possible for L not being of the order 2,3,5. If we 
know that far fewer than L/2 cells are active, then we 
can try to apply compressed sensing algorithms such as 
Basis Pursuit or Orthogonal Matching Pursuit to recover 
x from y = G(c)x. See 0, 0, iflOl , and 0 for 
descriptions of the recovery algorithms. 

In light of Theorem 12.81 we can extend Theorem 12.251 
in a different direction and obtain a large class of opera¬ 
tors that can be identified via regular operator sampling 
without knowledge of the support set. This class is larger 
than the class of area <1/2 considered in 0, 0. 

Theorem 2.26: Let T > 0, L E N and c € C L be cho¬ 
sen so that G(c) has full Spark, and let g = c n d n T- 

For 0 < A < 1, define the operator class FLt,l{ A) to 











13 


be the collection of operators H in OPW 2 (U 2 ) such 
that for some fixed fundamental domain R of the lattice 

(TL)7L x (1/T)Z , supp r\u = Sjj Q R and 

y! Xs H +(kT,e/(TL)) < A L a.e. (22) 

k,t 

Then the following are equivalent. 

(i) A < 1/2 + 1/(2 L). 

(ii) For all Hi, H 2 € Ht,l{ A), H\g = H 2 g implies 
Hi = H 2 

(iii) TLt,l{ A) is identifiable with identifier g in the 
sense of © 

Boelcskei and Heckel (0, 0]) have shown that, for 
operator classes like those in Theorem 12.251 if only L — 
1 cells are active, these can be determined, and hence 
the operator class can be identified without knowing the 
spreading support. Their analysis and derived recovery 
algorithms rely on the fact that by varying (f, v) you 
obtain a family of equations y(i, v) = G(c)x(i, v) where 
the vectors x(i, v) have identical sparsity structure. This 
allows for the recovery of almost every operator t lflOl . 
It5l l in the given class. 

In Theorem 12.81 we give up joint sparsity, i.e., the 
sparsity structure of x(i, v) varies with (i, v). A compro¬ 
mise based on the characterization found in Lemma [2.1 II 
that guarantees joint sparsity and allows us to use 
Theorem 12.81 is given by the following generalization 
of Theorem 3 in O (cf. Theorem 3 in 0). Note that 
the additional parameter K can be chosen independently 
of T and L, i.e., choosing K large does not increase the 
sampling rate, nor the size of the compressive sensing 
problem, i.e., of the matrix G(c). 

Theorem 2.27: For T > 0, and L £ N, let c £ C L be 
chosen so that G(c) has full Spark. Given K £ N, define 
the operator class TLt,l,k A OPW (M 2 ) by H £ TLt,l,k 
if and only if supp rjn = Sh C [0, LT] x [0,1/T] 
satisfies 

( a ) Xs H +(kT,e/(TL)) < L — 1 a.e., and 

k,l 

(b) each set Aj in the partition of [0, T] x [0,1 /(LT)] 
given in Lemma 12.1 II can be written as a union 
of sets of the form [0 ,T/K\ x [0,1 /(KLT)\ + 
(. qT/K,m/(KLT )), 0 <q,m< K. 


III. Proofs of Theorems 
A. Proof of Theorem 17.61 

Since S is closed, each f-section St of S is closed and, 
hence, measurable. Therefore, xs(7, •) is a nonnegative 
measurable function and J R x.s(t, v) do £ [0, oo] is well 
defined for all t £ R. It suffices to show the result for 
Ax> = || / r Xs('j u ) tfo'H finite, the infinite case then 
follows from this. 

Assume that A is a set of sampling with D(A) < 

^ Aqq. 

Then, we can choose a set P with positive measure 
and f R xs(t i v) do > cioo for all t € P. Assume without 
loss of generality P C [0,1]. For any e, there exist 
mt £ PW(St) with ||mt||i,2 = 1 and ||r7T,t| a ||^2 < e, 
t € P. Define kh{x, y) = m x - y (y ) for x — y € P, and 
0 otherwise. Then hn{x , t) = kh{x, x — t) = mt(x — t) 
and rjnit, v) = fntiy) for t £ P, and 0 otherwise, 
so H £ OPW 2 (S). Observe that [[cth/Il 2 = \/\P\- 
Note that it is easily seen that if XagA c a(>a identifies 
OPW 2 (S), then (c\) is bounded. Also, by hypothesis, 
there exists K £ N which bounds the cardinality of 
A 0 [x, x + 1] above for all x £ M. We compute 

\\ H T, c Mv= [ \^2c X KH(x,X)\ 2 dx 

AeA ^ AeA 

c\m x -\(X)\ 2 dx 

AeA 

<II(ca)II?. 

<II(ca)||| 

= ||(ca)||^ 

= \\(cx)\\l 
< \\(cx)\\l 


B. Proof of Equation ([9]) 

Definition 3.1: The non-normalized Zak Transform is 
defined for / £ <S(R), and a > 0 by 


^2m x - X (X)\ 2 dx 


AeA 


K y~] f \m x -x(X)\ 2 dx 
AeA ^ 

/*A+1 

K / \m x -x(X)\ 2 dx 

AeA^ 

K [ V] \m t (X)\ 2 dt 
Jo AeA 

K f 1 e 2 dt = || (ca) ||fo= K e 2 



Then almost every operator in H £ Ht,l,k can be 
identified by regular operator sampling. 

Note that alternatively to choosing K, we could at¬ 
tempt to introduce joint sparsity by assuming that, for 
example, r/n is smooth. 


Z a f(fo) = Y J f(t-an)e 2 ~. 

n£Z 

Z a f(t , u) satisfies the quasi-periodicity relations 

Zaf (t + a,u) = e 2niai ' Z a f (f, v) 
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and 

Z a f(t, v + 1/a) = Z a f(t, u). 

y/a Z a can be extended to a unitary operator from L 2 (R) 
onto L 2 ([0, a] x [0,1/a]). 

The following Lemma connects the output Hg(x) 
where g is a delta-train, to the spreading function 

Lemma 3.2: Let a > 0 be given and let g = E n &na- 
Then for all (t, v) £ R 2 , 

{Z a O H)g(t, u) 

= a^ 1 Y^ Y^ tjh (f + ak, u + m/a) e ^ 2muka ^ 

k m 

where yjj is the spreading function of the operator H. 

Proof: It can be verified by direct calculation that 
if 9 = E„^no then (Hg,f) = {r) H ,Z a f) for all 
/ £ <S(R) where the bracket on the left is the L 2 inner 
product on R and that on the right the If inner product 
on the rectangle [0, a]x[0,1/a]. Periodizing the integral 
defining the L 2 inner product on the left gives 

(' m,Z a f)= / Y, y: + ka,is + m/a) 

Jo k m 

e"2 niuka Z a f(t,v)dtdv. 

Since this holds for every / £ <S(R), the result follows. 


Lemma 3.3: Let T, LI > 0 be given such that TLl = 
1/L for some L £ N, let (c n ) be a period-L sequence, 
and define g = E n c n d n T- Then for (t, v) £ R x R, 

{Zi/n ° H)g{t,u) 

L—l 

= «E-.EE r]H{t + k/Ll + qT, v + mLl) 

g =0 k m 

27 ri(u-\-mQ)qT ^—2nii'k/Q 

(23) 

Proof: Note first that letting j = nL — q, 0 < q < 
L — 1, n £ Z, 

L -1 

9 = Y ( f S nT = YY CnL-q5nLT-qT 

q= 0 nEZ 

= Y c -Q T -q/LCl ( Y Sn f } 

q =0 ^ nGZ 

For a £ R, the spreading function of H o T a is r//y (t — 
a, v) e 2niva and hence 


{Zi/u°H)( Y c j 5 "-Tj{t, 


L -1 


Y c -q( z i/n oHo T -q/m) f Y Sn / n ) v )- 

q —0 '"nEZ 


Lemma [T2l yields the result. ■ 

Changing summation indices in (1231) by m = nL + £, 
0 < £ < L — 1, n € Z, yields the following lemma. 

Lemma 3.4: Let T, > 0 be given such that Til = 
1/L for some L £ N, let (c n ) be a period-L sequence. 
Then with g = E n SnT, and for all (t,v) € R 2 , 

(Zyn o H)g{t,v) 

L—l L—l 

= nY c -q YY 7 ] h F '(* + qT,v + mLl) 

q =0 m =0 k 

27 xivqT ^—‘I'Kivmq!L (24) 


where y^ P {t, v) is the (1/Ll, 1/T)—quasiperiodization 
of yu defined below. 

Definition 3.5: Given a bivariate function f(t, u) 
and parameters T, LI > 0, define the (1/Ll, 1/T)— 
quasiperiodization of /, denoted f® p , by 

f QP (f v) = YYH t + p + ^/ T ) 

k t 

(25) 

whenever the sum is defined. Note that f® p (t,v + 
1/T) = f QP (t,v) and f® p (t + 1/Ll,v) = 

e 2niu/Ci fQP( ti v ) for all (t, v) € R 2 . 

Lemma 3.6: Suppose that supp(/) = S is contained 
in a fundamental domain of 1/HZ x 1/TZ. Then 

fit, v) = YY f QP i f - p - l / T ) 

k t 

X [0 ,1 /n] (* - m) X[0,i/T\(v - t/T) X s (t, v) 

(26) 


where if / € L 2 (R 2 ), the sum converges in L 2 and 
uniformly on compact sets. 

Proof: Under the given assumptions, the functions 
being summed in (1261 ) have pairwise disjoint supports. 
Since 151 < 1, the sum converges in L 2 if / £ L 2 (R 2 ). 
Moreover, since on each compact set, the sum is finite, 
we get uniform convergence on compact sets. 

To complete the proof, we show that (l26l ) holds point- 
wise. Since S is a fundamental domain, for(f, v) € S 
only the (. k, t) = (0,0) term survives in (l25l ). Hence, for 
all (t, u), 

f QP (t,u)X s (t,u) = f(t,u). 


By direct calculation, 

f Qp (t,u) = YY f QP (t-k/Ll,u-£/T) 

k l 

X [0 ,t/n](t - k/Ll) X [0il/T] (u - t/T) e 2 ” ik ^ 

for each (t, u) £ R 2 . ■ 

Lemma 3.7: Let T, LI > 0 be given such that Til = 
1/L for some L £ N, let (c n ) be a period-L sequence. 
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Then with g = Yh n Cn $ nT ’ (*> v ) £ anc * P = 

0,1,-.., L-l, 

e -^T P{Zi/noH)g{t + Tpv ) 

L—1 

= (2 E (T'? M rn c) p e - 2mvT( i ri% P (t + Tg, z. + flm). 

q, m= 0 

(27) 

Proof: By (l24l) . 

(^l/Q °H)g(t + pT,u) 

L—l L-l 

= ^ E c -9 E i]% p (t + (q + p)T, V + m(2) 

5=0 £=0 

^—2'KivqT g— 2'rrii'mq/L 

Making the change of index q q—p, rearranging terms 
and using the fact that LT =1/(2 yields 

(•^t/n o H)g(t + Tp, v) 

L-l L-l 

—2'Kim(q—p)/L 

q= 0 m= 0 

7/g P (f + gT, v + m(2) e - 2 ™( q - p)T . 

Since (T q M m c) p = c p - q e ^^(p~<i )/ L , the result fol¬ 
lows. ■ 

Letting 

Z Hg (f v) p = (Z 1/n O JT) 9 (f + pT, v ) e- 2 ™ pT (28) 

and 



P iL (/, ^)(g,m) 


we have that 


(2 7?g P (f+^T, ly+mfl) e" 2 ™ 1 "^ e - 2 ™l m /L 

(29) 


L-l 

ZiHgift l y )p /* ' ^{^)p,(q,m) Vh (ti ^)(g,m) 

q,m=0 

which is ©. 


C. Proof of Theorem 12.61 

We first recall and outline the proof of Theorem 12.41 

Given any square submatrix of G(c), call it M, 
det(M) is a homogeneous polynomial of degree L in 
the L variables co, ci, ..., cl- i- In order to show that 
this polynomial does not vanish identically, it suffices to 
show that there is at least one monomial in det (M) with 
a nonzero coefficient. 

Such a monomial, pm, is defined recursively as fol¬ 
lows. If M is 1 x 1, then det(M) is a multiple of a 
single variable Cj and we define pm = Cj . If M is dx d, 
let Cj be the variable of lowest index appearing in M. 
Choose any entry of M in which Cj appears, eliminate 


from M the row and column containing that entry, and 
call the remaining (d — 1) x (d — 1) matrix M’. Define 

Pm = Cj pm'- 

The remainder of the proof consists of showing that 
the coefficient of pm is nonzero. In fact, it is a product of 
minors of Wl which, since L is prime, never vanish due 
to a classical result known as Chebotarev’s Theorem. 

Proof: (Theorem 12.61 ) Let k < L and choose 
k columns of G(c). Applying the algorithm described 
above to the non-square L x k matrix Mq formed by 
those columns, we can identify a monomial pm 0 - The 
key observation is that at each step in the algorithm, 
a variable Cj appears for which 0 < j < k. Once k 
rows of Mo have been eliminated, define M to be the 
k x k submatrix of G(c) consisting of those rows and 
the columns of G(c) chosen originally. 

Since the polynomial det(M) is not identically zero, 
and since at least one nonvanishing monomial of det(M) 
has only variables c } for 0 < j < k appearing, there is 
a c € C L , with supp(c) C {0, ..., k — 1} such that the 
columns of G(c) are linearly independent. 

Since the exceptional set of such c is the zero set of a 
polynomial in k variables, its complement is dense and 
open in C fc x {0}. Hence the (finite) intersection of these 
sets over all choices of k columns of G(c) is also dense 
and open in C fc x {0}. ■ 


D. Proof of Theorem [27 

Proof: Note first that by (l28l) and 

L 

|Z Hg(t, o) p \ 2 dt dv 


E 

P=0 7 7[0,T]x[0,fi] 


[0,T]x[0,LQ] 


\Zi/ n Hg(t,o)\ 2 dtdi' = \\Hg\\ 2 L2 


and 


E \rj H (t,v) q)m \ 2 dtdv 

q, m =oMo,T]x[o,n] 

= (2||r/EllL 2 ([o,LT]x[o,Ln])- 


(i)=Kiii). If (fill fails, then there exist integers qo and 
mo with S' = S (1 S+(moLT,noLTl) is a set of positive 
measure. This implies that there exists an operator H € 
OPW 2 {S) with spreading function rj € L 2 (R) \ {0} 
and r/2 p = 0. Indeed, as S' C S+(moLT, noLTi), we 
have S', S" = S'—(moLT,noLTl) C S and g(t,u) = 
XS'(t, v) — xs"{f^)c 2 ' Kwrn ° l/ 0 but rp p = 0. Then 
ZjL 9 = 0 which is equivalent to Hg = 0, showing that 
(i) fails. 
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Assume now that CED holds, and, without loss 
of generality, slippy C [0 ,LT] x [0,LQ], so that 

V® P \[o,LT]x[o,Ln] = V- 

If (IT 2 T) fails,Then there exists a set of positive measure 
A with Xs+(kT/n)(t, v)>L + l, ( t , v) € A. 

Hence, there exists A C A of positive measure and 
a fixed collection of L + 1 rectangles R /, ^ indexed 
by A out of the L 2 rectangles of size [0,T]x[0,H] 
tiling [0, LT] x [0, Ltt] with 22(k/)eAXs+(kT,m)(t,v] > 
L + 1, (t, v) € A. Since G(c)|a has L + 1 linearly 
dependent columns, we can choose a nontrivial vector x 
supported on A with with 0 = G(c)x, and, this allows us 
similarly to above to define a function t]h f 0 supported 
on A C S with G(c)rj H (t,u) = 0. As before, we 
conclude that Hg = 0 while H ^ 0. 

(iii)=Hii). Following the arguments above, 

A\\H\\hs < II-H^IIl 2 < B\\H\\ hs 

with A is the minimum over all singular values of Lx L 
sub-matrices of G(c) and B is the maximum over all 
singular values of L x L sub-matrices of G(c). 

(ii)=Md). Obvious. ■ 


E. Proof of Theorem \2.13\ 

Proof: Let L € N be as described. We will show 
that S meets at most aL rectangles R q ^ m , with T = y/~L. 
To this end, note that a Jordan curve Ci with length 
Ui £ ((fc ? ; — l)/\/X, ki/y/L), ki £ N, touches at most 4 ki 
rectangles R q: m , in fact, this bound is rather pessimistic 
and only sharp for k n = \. Note that 

N 

Vlu > vi 

i —1 

N N 

>VLj2(ki~l)/VL=(J2 k i)~ N , 

1=1 1=1 

and, hence, the number of rectangles B{dS ) needed to 
cover the boundary dS of S satisfies 

N N 

B(dS ) < B (Ci) <Y 4k ^ ^(VlU + N ). 

i= 1 i=l 

We conclude that the ”fat” boundai'y, that is, the 1 /\[L x 
1 /VT rectification of the boundary has area bounded 
above by 

4 (yfLU + N)/(yfLf = A(U/yfL + N/L) < e. 

It follows immediately, that at most oL sets R, hm are 
needed to cover S. ■ 


F. Proof of Lemma \2.11\ 

Proof: Note first that with S° given by (fl3l) . 

Yj XS+{kT//(TL)) = Yj XS°+{qT,m/{TL)) 
k,£ 0 <q,m<L 

so that (fl 2 l) is equivalent to 


SUP Y XS°+{qT,m/(TL)){t,v) < L. 

(t,y)e[0,T]x[0,t/(TL)] Q< q m<L 


Assume that (fl2l) holds. Then for each (t , u) € [0, T] x 
[0,1/(TL)] there is a unique 0 < n < L and |T| = n 
such that 


Yj XS°+(qT,m/{TL )) v) 

0 <q,m<L 

= Yj XS°+{qT,m/{TL)){t,v) = n. (30) 
(g,m)er 

For each such n and T, define the set 


A n ,T = {(t,v) £ [0, T] x [0,1 /(TL)\: ([HDD holds}. 


This collection of sets forms the desired partition of 
[0,T] x [0,1/(TL)]. It is clear that the sets Sj defined 
in 03} satisfy the required conditions. 

For the other implication, if Aj + ( kT, f/(LT)), 0 < 
k,£ < L, meets S° at most L times, then 


sup 

0,v)e[0,T]x[0,l/(TL)] 

sup 

(t,^)e[0,T]x[0,l/(TL)] 


Yj XS°+(qT,m/(TL)) v ) 
0 <q,m<L 


sup , E XAj+(qT /(TL))(t,o) < L. 

1 - J - JV 0 <q,m<L 


G. Proof of Theorem \2.15\ 

Proof: Suppose first that (foj^o) = (0,0), and that 
c is chosen so that G(c ) has full spark. By the support 
assumption on S, (1271) implies that for 0 < p < L — 1, 
© takes the form 


L-l 

Z Hg{t , v)p = E G M», (S ,m) VH(t,v)(q ,m) • 

q,m =0 

Let [l>( q , m/) . p \ be a left-inverse of the L x |r| matrix 

[ G ( c )p,(9,m)]o<p<L,(9,m)er- That is, for every (q,m), 

(■q',m') € r, 


L—1 

Y^j b(q,m),p G(c) p ^ q r m /' ) 

p =0 


1 

n 


e 2n iqm/L^_^ 


5 


m—m' • 


( 31 ) 
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Again by support considerations on S, i/n E 
OPW 2 (S ) satisfies 

?7^ P (f,i/)X[o,i/n](t)X[o,i/T](^) 

= XI l/ ) x [o,r](i - q r T’)X[ 0 j n]- run), 

(q,m)e r 

and for each (q, m) € T, 

V h P (t + qT,v + mn)X [OT] (t)X [on] (i^) 


L -1 


0 2mv(q-p)T 


= X e 2 ” fci//n E ^ N p 2iri(v-£/T-mtt)(q-p)T 


= 2^ e 

k,e 

—2'Kivk^l 


k,l 

(Z 1/Sl O H)g(t - k/n -(q- p)T, u - t/T) 

[ x [ 0 ,T](t - k/£l - qT)X[ 0>Q] (u - t/T - mn)X s (t, u)\ 

L—l 

E 2mku/n t 2m(L'-mCl)(q-p)T 

e w W, m ),P e 

(g,m)er,p =0 

{Zi/n ° #)0(* - (? - p)r, v) 

[ x [ 0 ,T](i - k/Q - qT)X [o n] (v - t/T - mQ)X s (t , z/)] 
L—1 

X k(g,m),p e 

(g,m)er,p =0 

(^l/n ° - (9 - p)?\ z/) 

^2 x [0,T](t ~ k/n - qT) 
k,i 


0 2iri(u—mQ)(q—p)T 


x [o,n](v ~ V-/T - mti) X s (t, v) 


Defining 


s {q ,m) = s n (\J R q , m + (k/n, i/t)\ , 

'kl ' 


it follows that S = (J ( 9 m)er that the union is 

disjoint, and that 

X S iq , m) (t, ,) = X X [0 ,T] (t ~ k/n - qT) 


kl 


x [o&]{v - t/T - mtt) X s (t, v). 


Therefore, 


— X k(g,m),p x [o,r]( i ) x [o,fi]( 1 ') 

p=0 

{Zi/n o H)g(t + pT, v). 

Therefore, by the quasiperiodicity of the Zak transform, 

V% P (t, iy) x [o,i/n] (t)X[ o,t/T] (v) 

L -1 

= X X X [ 0 ,T](t - qT)X [0j n] (z/ - mQ) 

(q,m)er p =o 

e 2 ^-mmq-p)T (Zi/n 0 ^ )5( t - ( g _ p )T, v). 

Applying d26l ), 

m(t, v) = Vh P fa z/ ) x s(^ 

= X g% P (t-k/n,v- lIT) 

k,t 

IX,i/f 2 ](f - fc/O)X[ 0 ,i/r] -t/T)Xs(t,v)\ 

L -1 

'), , e 2 
q,m),p ° 

(g,m)€r,p =0 


0 —2 , Ki(y—m£l)(q—p)T 


L—l 

~ X b (<h m ),P e 

(g,m)el\p=0 

X w - n/n - (9 - p)T) e 2 ^ n Xs ( ,. m) (i,«/) 

nEZ 

rL-1 

= X X X e 

(q,m)er Lp= 0 nez 


27 ri(is-mQ) (q—p)T 


Hg(t — nLT — (q — p)T) e 


2niunLT 


X S(n,m.) (t, v) 


Extending b^ m ^ p to have period L in p, it follows that 


L—l 


—27 Ti(i/—mQ,)(q—{p—nL))T 


- E EEw.^ 

(q,m)er p=0 nSZ 

H 9(t - (9 - (P - « L )) T ) X S(q, m ) (t, v) 

= X X 6 fe-)^ e_27ri( " _mn)(9_fe)T 

(g,m)er fc 

Hg{t - (q - k)T ) X S(g m) (f, v). 

Finally, writing 

*„«) = /,<*, 

yields (fl 6 l) with (foj^o) = ( 0 , 0 ). 

To complete the proof, note that for almost every t, the 
set, {z^: (f, z/) E S 1 ^ ^} is contained in a fundamental 
domain of the lattice TZ of M. This implies that the 
measure of each such section is no more than 1/T, and in 
particular that for almost every t, Xg^ qm) (t, •) € L 2/ 
Therefore, by Plancherel’s Formula, 


I *h(q,m) ili s) | 2 dt ds 


„ 2niiys 


x S(„, m) (t, V) dv 


5(q,m) 
12 


ds dt 


\ x S (q , m) (t,v)\ 2 dudt = \S (qjtn) \ 2 < oo 

and for almost every (t, s ), 

l $ (g,m)(^s)l < J X S( q , m ) (t, 12 ) dv < 1 /T. 
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Hence <I ) ( r/ rn ) E L 2 n L°°(R 2 ). Convergence of the re¬ 
construction sum in L 2 (R 2 ) follows from the observation 
that Hg E L 2 (R) (see Lemma 13.21 ) and basic properties 
of the Zak Transform (see e.g., fill . Section 8.2). 

If (to,i/o) / (0,0), we formally compute 


combining (1331) with (|32T ) yields 

h(x, t) 

L -1 

EE b (q,m),k (M- Vo T- to H)g(t - to-{q - k)T) 

k j =0 

^—2-Kim.(q—k)/L ^—2ni(x—t—to+(q—k)T)v 0 


H = 


IS-(to,Vo) 
= Tt 0 M Uo H, 


M u T t dtdv 

Vh( t + t 0 ,v + v 0 ) T t+to M v+Vo dt dv 


$(q,m)(t, (x - t - to) + (q - k)T) 

_ e 2m(x+t—to)v 0 
k (q,m)eT 

^—2'Kim(q—k)/L ^-27ri(x—t—to-\-(q—k)T)uo 


where rjg(t,v) = e ~ 2mtu ° r)n{t + to,v + vq). Taking 
inverse Fourier transforms v —>• x on both sides, we 
obtain hg(t,x ) = hn(t + to, x) e~ 2mu ° x which 

is 


h H (t, x ) = e ^+t-to>o hs ( t _ tQj x y (32) 


®(q,m)(t, (x - t - to) + (q - k)T) 

_ p 2wi(t+to)vo 

E E b ( q ,m),k Hg(t — (q — k)T) e ~ 2 ^ im (<i- k )/ L 

k (q,m)er 

$(q,m)(t, (x - (t + to) + {q~ k)T ). 


With S = 5— (to, k'o), we can apply (fl6l) with (to, v o) = 
(0,0) to reconstruct h^ from Hg with the same g = 
X) c n 5 nT , that is, 

t ) 

= E E b Mk Hg(t-(q-k)T)e- 2 ^- k V L 

k (q,m)e T 

$(q,m)(t,(x-t) + (q-k)T). (33) 

where 

$(q,m) (t, s) = J f?™ S (t, V) dv 

and 

S(q,m) = sn \J (R q , m + (k/n,e/T)). 

k/GZ 

Observing that S( qjm ) = <S( 9im ) + (to,^o), we obtain 
®(q,m)(t, «) = e 27r “ V ° $(,, m )(t - t 0 , s). 


H. Outline of Proof of Theorem 12. 1 71 

The proof follows that of Theorem 12.151 once we es¬ 
tablish that we can replace the shaip cut-offs, Xr 0 j T 1 and 
X[ 0 i n] by smooth ones. Since S is compact and |S| < 1, 
for 5 > 0 sufficiently small, the set S$ = 5+[— 6, <5 ] 2 also 
satisfies | S$ \ < 1. Since Theorem 12.151 allows us to shift 
the region, and since IS) < 1 , we can assume without 
loss of generality that there exist T > 0 and L E N 
such that S C (0 ,TL) x (0,1/T) and that S$ has a 
(T, L)-rectification. Since S C U( 9 im ) e r R q ,m = R, it is 
sufficient to prove the theorem with OPW 2 (S) replaced 
by OPW 2 (R). 

By Lemma [3771 given H E OPW 2 (R) with spreading 
function ’quit, v), and given any weighted delta train of 
the form g = ff ri c„ S n r where c = (c„) is a period- 
L sequence, (l27l) holds with t)® P replaced by tjh for 
all ( t,v ) in an e-neighborhood of [0, T\ x [0, fi], R^ 0 = 
([—e/2,T + e/2] x [—e/2, fi + e/2]. 

Let r, ip E <S(R) satisfy 

suppr C [—e/2, T + e/2], (34) 

suppy? C [—e/2, H + e/2], 

so that suppr(f)^(i/) C _Rq 0 , and 

El r (t + kT) = 1 = El p{ v + nVt), (35) 

kGZ nGZ 

for all (t,v) E R 2 . For e < 5, it is not hard to show that 
if R q , m % R then 

r]H(t, v)r(t — qT)(p(v — mtt) = 0. (36) 
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Therefore, 

m(t,v)= Y ri% P {t,u) r {t-qT)(p{u-mLt). 
(g,m)er 

Following the proof of Theorem l2.15l with r(t) replacing 
x [o,T](t) and (p{u) replacing X [0 ,n]O)> 

m (t,u)= Y Y b M,k^ i{v - ma){q - k)T 

(q,m)G r k 

Hg(t - (q - k)T) R^ m ) (t, u) 

where 

P(q,m) (^> 1/1 ) 

= r(t — k/Ll — qT)(p{u — i/T — mil) X^(f, i/) 

k,t 

= r(f — gT) <^(z/ — mfl). 

Finally, 
h(x,t ) 

= E Y b M,ke 27 rimiq - k)/L Hg(t-(q-k)T) 

(q,m)e r fc 

(x-t) + (q- k)T) 

where here 

$(q,m) (t, S) = J e 2mUS R(q, m ) (t, u) du 

= r(t — qT) e 2nisrn ^p(s). 

Plugging this into (fl6l ) gives the result. 

I. Lattice tilings and proof of Theorem \2.18\ 

In this section we will prove Theorem 12.181 but also 
derive results where the tiling of S is defined by arbitrary 
full rank lattices in M 2 . The reconstruction formulas use 
results from representation theory; these carry over to the 
higher dimensional setting if the lattice is symplectic. 

Proof: As before, we assume that SCI 2 satisfies 
| S'| < 1. Suppose that for some A = ( 021022 ) with 
det A = 1/L, S is contained in a fundamental domain 
of the lattice LAI?. The lattice LAI? is the so-called ad¬ 
joint lattice A° of A. Indeed, A° = (1 (\[L A)° = 

yfLyfLA = LA (see fill for details). We shall assume 
without loss of generality that an f 0. Otherwise, we 
could replace the first column with the second and the 
second with the negative of the first, leading to a different 
parametrization of the same lattice. Further assume that 
there exist to,UQ, and T C Z 2 , |T| < L such that with 
Pq,m = ^4([0, l] 2 + (t 0 ,is 0 ) + (q, m) T ), q, m € Z, 

LA[ 0 ,l] 2 n ( [j S + LA{k,l) T 

X k,l£Z 

(37) 



As before, we will set 

$( 9 ,m) {t, s) = J e 2mus Xs (q , m) (t, v) du 

where 


S(q,m)=Sn 1J (P q , m + LA(k,£) T ). 

k/GZ 

We will derive reconstruction formulas and show that 
if 012/011 is rational, then OPW 2 (S) can be identified 
with a weighted delta train and if 021 an is rational as 
well, then the coefficient sequence (c n ) of that delta 
train is periodic and we are in the framework of regular 
operator sampling. 

We shall assign to each operator H € OPW 2 {S) 
an operator in H € OPW 2 (L~ l ?A~ l S) and then 
apply the reconstruction formula in Theorem 12.151 to 
reconstruct h = hg of H € OPW 2 (L~ 1 ?A~ 1 S). From 
this, we will construct h = hu and therefore H. 

The result is based on the existence of the operators 
p{y/~LA) that appear in the following computation. The 
existence follows from the representation theory of the 
Weyl-Heisenberg group and is discussed in this setting 
in El, El. Let p{t,u) = e™ tv T t M v , rj*{t,v) = 
e - mtu i](t,u), and B = \j~LA. Then 


H = 


rj{t, u) T f M v dtdu 
pit, 1 ;)e~ nitv e nitv T t M v dt du 
u) p(t, u) dt du 
rft 1 {B{t,u)) p{B{t,u)) dtdu 

B~HS ) 

r]#(B(t,u)) p{B)p{t,u) p(B)* dtdu 


= p{B) JJ p#iB(t,u)) p(t,u ) dtdu p(B)* 

= p{B ) H piB)* , 


with rj#{t,u) = q#{B(t,u)). Setting Qi{t,u) = t and 
Q 2 (t, u) = u we have 

77 (t,u) = e ni{ - tv -Q' B{ ~ t ' v kQ*B^)) r] (B{t,u)). 

Moreover, observe that S = B~ 1 S satisfies the 
hypothesis of Theorem 12.151 with T = Q = 1 /sfL. 
We have therefore with an L periodic sequence (c n ), 
9 = I?Cnd ny /Z’ and B ~ l = (-621 "&?”) the reconstruc¬ 
tion formulas 


h(x,t) = e 2 ^ t+t °>°Y E Krn),k 

k ( q,m)GT 

Hg{t -(q- k)/VL)e- 2nim ^~ k ^ L 

5 ( q ,m) (t+t 0 , x - (t+t 0 ) + (q- k)/y/L )], 
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k (<j,m)er 

Hg(t -{q- k)/VL)e~ 2 ™ m ( q - k V L 

(t+to, ^) e 2-(*+*o-( 9 -fe)/VI)-] 

_ e m(tv-QiB(t,v)-Q 2 B(t,v)) ^j,)) 


r)(t,v) = lB 1 (*’")+ t °) I/ ° 

e -iri(Q 1 B- 1 (t,v) Q 2 B- 1 (t,u)-tu)'^~^ ^ 

fc (g,m)sr 

Hg{Q\B~ l {t, v)- (q- k)/VL) 
e-™^ q - k VL XSi ((t,v)+ B (to, 0 )) 

e 2iri(Qi_B - 1 (t,v)+t 0 -(q-k) / y/T)Q 2 B~ 1 

_ e 27 r *((6 22 t-fe i 2 ^)^ o +i°(6 ii ^-6 2i t)+to^o) 

e ni((b 22 t-b 12 v)(b x xt'-b 21 t)-tv)'^^ ^ [fe(g m) fc 

fc (g,m)er 

Hg{{b 22 t - b 12 v) - (q- k)/VL) 

e - 2 «im( q -k)/L Xs ^ ^ v)+{blltQi 6 21 t o) ) 
g— 2wi(q— fe)/vT)(6iii/-62it)j 

_ g27rt ( [( fl 22^ — a 12 ( a ll ^— a 21^)] VT+to^o) 

^■Ki{L{a 22 t-a 12 u)-{a 1 - i _v—a 21 t)—tv) E E [W 

k (q,m)er 

Hg((a 22 t - ai 2 v)VL - (q — k)/\fL) 
e - 2 *im(q-k)IL Xs ^ + „ + vWt 0 ) 


= — 2 m{q— fc)(an!/—a 2 i 


* } ]- 


(38) 


which leads to 


_ g— 2m(\/Lt 0 i/o+L/2)a 22 a 21 t 2 ^2 -kUqVo 

E E [&(g,m),fc-^5(\ / I(a 2 2t - (? - fc)/L)) 

k (qf,m)Gr 
g—27 rim(q—k) / L 


$( 9 ,m)(i + VLant 0 t, x + fo^oair/k - (9 - &) a n) 

g— 2 - 7 riv / ia 2 it 0 ( 3 :+toi / oCiiiV / i—( 9 —fc)an) e 2m(q—k)a 2 it j 


(40) 


and, if fo = 0 , 

/i(x, t) = e - ,riLa22a2lt2 

E E [ b ( q ,m),k Hg{\fL(a 22 t - (q - k)/L)) 

k (q,m)e T 

g— 27 rim(q—k)/L ^ 27 ri(q—k)a 2 it 

x- {q- k)a u )] (41) 


By construction, we have Hg = n(B)*Hg(B)g with 
g = Hence, we can replace H in (l38l) by 

g{B)*H and g by g where g = /j,(B)g. In the following, 
we will give explicit representation of g(B) and examine 
g = fiiB jg. Note that the given reconstruction formulas 
hold true for any tempered distribution g = g(B)g, 
but we are mainly interested in the case that g(B)g is 
discretely supported, or, better, g = n(B)g = ^ CnbnT 
for some T > 0 and a periodic sequence c = (c n ). In 
applications, this would allow us to use any hardware de¬ 
veloped to excite an operator described in Theorem l2.15l 
Recall that B = \[LA , so det B = 1 and we assume 
611 / 0. We have 


Taking inverse Fourier transforms u —>• x on both 
sides gives us a formula for h, but as the right hand 
side contains the product of three functions in u, the 
resulting formula for h does not give much insight in 
general. If a \ 2 = 0 though, the above simplifies (using 
<211022 = 1 /L) tO 


V(t, v) = '^2 Y [ b (Q^),k Hg(a 22 VLt - (q - k)/\fl) 


( On 0 12 \ 

V 621 i >22 / 

= (a il ?)(?-o 1 )(-^ 1 2?)(- 0 iJ)( b o 1 i/° bll ) (42) 

Using notation from ifllll . we have 

/ 21 (a) = /r(i J) :/■->• e mabr f, 

r = n{°ih) 

/22(a) = / 2 (o i/a) : / a _1/2 /( • / a), 


k (q,m)£ r 


hence, 


e - 2 nim(q-k)IL ( * + „ + VZa 21 t 0 ) 

g 27 ri(l 0 v 0 aiiVT—( 5 —fc)an)i/ e ~2ni(\/Ltoi'o+L/2)a 22 a 21 t 2 
^2-n-itoUo e 2m(q—k)a 21 t 



m({> : 


bn b\ 


) = gi{b 2 i/bn) B*gi(—bnbi 2 ) Tg 2 (bu) 
= / 2 i(a 2 i/an) .FVi(-Laiiai 2 ) -E/^(V^an). 
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This leads to 

v(B)g = b“ III ) X] Cn ^n/VL 

= ^i(a 2 i/a n ) J 7 *Hi(-Laua 12 ) Bo 

^(VTan) E Cn ^n/VL 


if a 2 ian is rational as well, then we are assured that the 
coefficient sequence (c n ) has period 

L" = LCM{q[a 21 an(L/L') 2 /2], L'} 

= LCM{q[a 21 an(L/q[ai 2 /(2Laii)]) 2 /2], 
q[ai 2 /( 2 La n )], L}, 


= (VToii) 1/2 /Ui(a 2 i/a n ) T*fii(-Lanai2)o 
B ^ ^ Cn^nau 


= (\/La n ) 1/2 qi(a 2 i/an)o 

Lanai 2 ) E c m^m/(Lan) 


= (vTa n ) 1/2 /Ui(a 2 i/a n )o 

-2nim 2 a 12 /(2La 11 ) c 

J / j c, m/(Laii) 


that is, we are in the framework of regular operator 
sampling. 

Let us consider the special case that oi 2 /(2an) is an 
integer (for example, if a \ 2 = 0 as in Theorem 12.181) . 
then q[ai 2 /(2Lan)] £ {1, L}, so L' = L and L" = 
LCM{q[a 2 ian/2], L}. If in addition La 2 ian/2 is an 
integer, then q[a 2 ian/2] € {1 ,L} and L" = L. 

To complete the proof of Theorem 12.181 observe first 
that L = U , and indeed (c n ) = ( c ' n ). Consequently 

g = fi(B)g = Y J c n e^ in2a ^5 nail . 


where we have used the fact that the Fourier transform of 
a delta train of the form ^ ngZ c n^nT, where c = ( c n ) 
has period L is another delta train of the same form. 
Specifically, 

B ^ CnSnT = ^m/LT (43) 

nEZ mEZ 

where c denotes the Discrete Fourier Transform of c, that 
is 

C m = C k e-^n/L. 
k =0 

Equation (1431) is a simple consequence of the fact that 


Further, observe that 

, (fen b 12 \* 

FU 21 fe 22 ) 

= j u 2 (>/Iaii)*^/xi(-Loiiai 2 )*^/x 1 (a 2 i/oii)* 

= /i 2 (l/(A/Zan)) /ii(-a 2 i/an). 


Hence, if ai 2 = 0, then 

KklYfW 

= /r 2 (l/ (VLan)) ^i(-a 2 i/an)/(x) 

= (VZan) 1 / 2 e -™2i/an(VIa ll3: )= y(VZa n x) 
= (VZan ) 1/2 e-^-anx 2 f(VLa u x) 


Fj2 S nW = ^ 

n SZ m 

The sequence e ~ 2mm2a ^ /( 2 Lan) j s p er j 0 ( ji c j n m jf 
e -27nmai 2 /(2Lan) ff, a f j Sf jf ai 2 /an is rational. In the 
following, LCM refers to least common multiples of nat¬ 
ural numbers, and for a rational number a, q[a] denotes 
the smallest natural number q such that qa is an integer. 
With this notation, (c') m = c m e ^ 2mma ^/{ 2 La 11 ) f orms 
a sequence with period L’ = LCM{q[ai 2 /(2Lan)], L}. 
Once again employing (l43t . 

g,{B)g = (\/Ian ) _1/2 gi(a 2 i/a n ) T* ^2( c ’) m 5 m /( Lail) 
= (\/Ia n ) _1/2 gi(a 2 i/an) ^2c' n 5 nailL/L , 

= (VLan)- 1/2 E c n eWa2iail(i/L,)2/2 ^aiiL/L'. 

(44) 


and 


g(B)*Hg(VL(a 2 2 t ~(q~ k)/L )) 

= (v/Zan) 1 / 2 e -’ rii “2i a ii (VL(a 22 t.-(q-k)/L)) 2 

H 5 (\/ZanV / Z(a 22 f — (q — k)/L )) 

= (VTan) 1 / 2 e -’ ria 2i“i i(La 22 t-(q-k)) 2 

Hg(t - a u (q - k )) 

We conclude that 


= (/La n ) 1 / 2 e- 2 ^ (VIt 0 l/ 0 +L/ 2 )a 22 a 2 lt 2 e 2,rit0l/0 

(q,m),k & 


y ' ^ N y f> —ma 2 i_a 1 i(La 22 t—(q—k))' 2 


k (g,m)er 

iTs(t - an(q - A:)) e -^ m ^~ k )l L 


$( q ,m)(t + VLa u t 0 t, x + t 0 i 2 0 a u VL - (q - k)an) 


We conclude that g(B)g = Yh c n ^nT with T = 
anL/q[ai 2 /( 2 Lan)] if ai 2 /an is rational. Moreover, 


e —2iriy/ r La, 21 to(x+t 0 v 0 a 11 y/L—(q—k)a 11 ) e2 iTi(q—k)a 21 t^ 

( 45 ) 
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and, if to = 0, 


h(x,t ) = (\/Lan) 1//2 
b(q,m),k e 


g — TriL(l220,2lt 2 

-7ria 2 iaii(La 2 2t—(<?—fc)) 2 


fc (g,m)Gr 


ff^(t - On(g - fc)) 


(46) 


.]. Proof of assertion in Example 12.201 

The goal is to show that OPW 2 («S) where 5 is the 
region shown in Figure [5] cannot be identified by regular 
operator sampling for any T or L. We will show that 
the TZ x (1/LT)Z periodization of S' does not form an 
exact L-cover for any T or A, thus violating (fl2l) and 
Remark 12.91 2) 

Proof: Assume first that T is rational. We can 
assume without loss of generality that T = 1/A' for 
K E N, and hence that Q = l/LT = K/L is also 
rational. Indeed, if T = p/q and if for some L E N, 
g = Yh c n^nT ( c n with period L) identifies OPW 2 (S) 
then letting V = l/q, L' = pL, and d n = c n / p if p 
divides n and zero otherwise, then d n has period L' and 
g = dribnT'- Note that the set of discontinuities of 

the function Xg + t kT en)> (k,£) € Z 2 in the rectangle 
A=[ 0 ,T] x [ 0 , Q] must occur on line segments of slope 
// or ^ \ passing through R (that is, intersecting 

two edges of R). In order that Y X-s+ikT/Ct) be constant 
and hence continuous on R, each such segment must 
coincide with at least one (kT,£Sl )-shift of a different 
such segment. In particular, the segment of slope / 
containing (0,0) and intersecting one side of R must 
be met by some ( kT , ^fi)-shift of the segment joining 
(2, y/2) and (4,2\/2 + 1/2), which implies that this 
segment must contain a point of the form (kT,£Q). 
However, a simple calculation shows that since T and 
O are rational, this is impossible. 

Now assume that T, and hence also H = 1/(AT) 
is irrational. In this case, discontinuities of Xs+(kT,en)i 
in the rectangle R must lie on lines passing through R 
with slopes as above, or on a pair of line segments of 
those slopes terminating at their intersection point in the 
interior of R (see Figure [ 6 |b)). There are at least one 
and at most three shifts with discontinuities of the latter 
type. To see this, note that since T is irrational, neither 
( 2 , y/2) nor ( 2 , s/2 + 1 / 2 ) lies on a vertical grid line of 
the form t = rriT and that since f> is irrational at least 
one of these points does not lie on a horizontal grid line 
of the form a = nil. Similarly, (4, 2\/2 + l/2) cannot lie 
on a vertical grid line but may lie on a horizontal grid 


line. Then by considering the cases in which exactly 1, 
2, or 3 of these points do not lie on a horizontal grid 
line, it is clear that in order for all discontinuities to be 
resolved, (4, 2\/2 +1/2) must differ from either (2, s/2) 
or (2, \/2 + 1/2) by some (kT,£ fl), which is impossible 
since T is irrational. ■ 


K. Proof of Theorem 12.221 

Proof: Let L > N 2 be prime, let Q = 1 /(TL) and 
let R q) m , {hi m) € T be the rectangles in the (T, N)- 
rectification of S. Then letting 


R 'q',m’ = [0, T] x [0, H] + ( q'T,m'n ), 

(q > ,m!) E Z 2 , each rectangle R q: m is covered by a 
collection of rectangles R' q , m , satisfying 


E 

{(g',m'): R q , 




f \Rq,m\ 



Let T' be those ( q',m') € Z 2 such that R' q , m , has 
nonempty intersection with S°. Therefore, 


\Rq',m'\ 

S E E 

(q,m)e r {( q',m '): R', ,C -R q ,m/ 0 } 


I R 1 I 

I -^q' ,m' I 


Rq,m\ 


S E 

(g,m)e r 

„ , 2JV 

— \ R q,m\ H J— 


(q,m)er 

irl 2 N 

< — H-< 

~ N L ~ 


L 

|r| +2 

N 


<|S|(l + e). 


Consequently, |r'|/A < |5|(1 + e), and S° C 
U(q/ im /) g p/Aq, m , = R. By Theorem 12.61 we can choose 
c € C L such that ||c||o < |T'|, Spark(G(c )) = |r'| and 
c is supported on its first ||c||o indices. Since S° C R, 
any identifier of OPW 2 (R ) is also an identifier of 
OPW 2 (S). Since R consists of only |r'| rectangles, it 
follows that vector on the right side of ([fi]) has at most 
|r'| nonzero entries and hence is solvable as long as 
Spark(G(c )) = |r'|. From this it follows immediately 
that Yn c n/nT identifies OPW 2 {R) and 

**<S<|5|(l + e) . 
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L. Proof of Theorem 12.251 and Theorem 12.261 

Proof: (Theorem 12.251) By Theorem 12.131 we 
can choose L € N so that every operator in 
TL(A, B, U, N, e, 1/2) has the property that suppp 
touches at most L/2 sets of the form 

Rq, m = [0, 1/VI] X [0, 1/VI] + (g/VI,m/Vl), 

(47) 

q,m = —(L - l)/2, —(L - l)/2 + 1,..., (L - l)/2. 

Now, let {S m : m = 1,, (V)} be the collection 
of area 1 sets that are formed by exactly L subsets of 
the form R, hrn in (l47l) . Choosing c € C L so that G(c) 
is full spark, it follows that for each m, OPW(S m ) 
is identifiable with identifier XVez ° n VtCT an< ^ that 
constants C \, C 2 > 0 exist such that 


This is easily seen by considering the sets 

\J[R q ,m + (kTL,£/T)\nR 

k,i 

where for each 0 < q,m < L, 

R q ,m = ([0, T\ x [0, 1 /TL\) + (qT, m/TL). 

Then S\ and S 2 can be formed by choosing two disjoint 
collections of \(L + l)/2] such sets. Since 

Y, X(s 1 uS 2 )+{kT,e/(TL)) > (L + 1) a.e. 
k,t 

the same argument as in the proof of Theorem |2^8] allows 
us to define distinct operators Hi, H 2 G 'Hr.rfA) with 
supp t]h 1 C S i and supp ?/# 2 *= $2 such that (Hi — 
H 2 )g = 0. Hence (ii) fails to hold. ■ 


Ci\\H\\ hs < WH^K/VlWv < C2\\H\\hs, 

nGZ 


for all 

He J OPW 2 (S m ). 

The proof is complete by observing that for H \ . H 2 € 
TL(A, B, U, N, e, 1/2) (which is not a linear space), we 
have Hi — H 2 € OPW 2 (S m ) for some m, and, hence, 

- H2 \\hs < \\(Hi - H 2 )Y, c nS n/VI \\L> 

< C 2 \\Hi — H 2 \\hs, 


Hi, H 2 g R(A, B, U, N, e, 1/2). Clearly, this leads also 
to the weaker statement (Hi — H 2 ) n = 0 

implies Hi = H 2 . ■ 

Proof: (Theorem 12.261) The proof of this result 
follows the proof of Theorem 12.81 
(i)=Kiii) Note first that if A < 1/2 + 1/(2 L) then 
A L < (L + l)/2. Hence if Hi, H 2 G TLt,l(A), then 
with supp r) Hl -H 2 ^ suppr/#, U supp ?/# 2 = S, 


Y, XS+{kT,l/{TL)) 
k,t 


< 2A L < L + 1 


and since the left side of the inequality is an integer, 
Hi — H 2 G H-t,l( !)• Therefore, (fTTT) and (fl2l) hold, and 
by the same argument as in the proof of Theorem 12.81 
(iii) holds. 

(iii)=Hii) Obvious. 

(ii)=>(i) Suppose that A > l/2 + l/(2L). Then we can 
find disjoint sets Si, S 2 C R such that 


Y XSi+(kT,e/(TL)) 
k,e 


oo 


A L > (L + l)/2 


IV. Conclusion 

This paper contains results relevant to two questions 
on the identification and recovery of operators with 
bandlimited symbols from the response of the operator 
to a regular weighted delta train. Such operators model 
time-variant linear communication channels. When the 
identifier is a weighted delta train, we refer to this the 
identification as operator sampling and when the weight¬ 
ing sequence is periodic as regular operator sampling 
The procedure is a generalization of classical sampling 
results for bandlimited functions, and of the determi¬ 
nation of a time-invariant communication channel by 
measuring its response to a unit impulse. 

We obtain a simple condition on the set S that char¬ 
acterizes when OPW 2 (S) can be identified by regular 
operator sampling. The condition requires that S be 
contained in a fundamental domain of a rectangular 
lattice and that its periodization on a reciprocal lattice be 
bounded above by a constant depending on the lattice. In 
this case, |5| < 1 , and we obtain explicit reconstruction 
formulas for the operators in OPW 2 (S). We consider 
the case in which S is contained in a fundamental 
domain of a general symplectic lattice and give sufficient 
conditions on the lattice under which OPW 2 (S ) can 
be identified by regular operator sampling and obtain 
explicit reconstruction formulas in this case as well. We 
provide an example of a set S for which OPW 2 (S) can 
be identified by operator sampling but not by regular 
operator sampling. 

For these results it is required that the support set 
be known. We also obtain a result showing that, under 
mild geometric conditions, recovery is possible when the 
support set is unknown but has area smaller than 1/2 and 
we characterize all support sets for which identification is 
possible via regular operator sampling when the support 
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set has area < 1/2. This characterization allows us to 
define a large class of operators for which identification 
is possible when the spreading support is small. This 
class includes the class similarly characterized in 01. |[4l. 
It is shown in [31, 0 that this class can be identified 
without knowledge of the spreading support for areas 
less than one. Following the ideas given in |3j], 0, we 
define a larger class of operators with area less than 
one that can be similarly identified without knowing the 
spreading support. 

Finally, we give a necessary condition on the rate of 
sampling, that is, the average number of deltas in the 
identifying weighted delta train per unit time, required 
to identify an operator with bandlimited symbol. The 
necessary rate depends on the bandwidth of the spreading 
support. We give a sufficient condition on the sampling 
rate in terms of the area of the spreading region. As 
a consequence of this result, it is observed that if the 
area of the spreading support is small, then any operator 
in the class of operators having that spreading support 
can be identified by only a portion of its response to 
an appropriate identifier. The fraction of the response 
sufficient for identification is asymptotically proportional 
to the area of the spreading support. 
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